Universality classes in creep rupture 
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We study the creep response of solids to a constant external load in the framework of a novel fiber 
bundle model introduced. Analytical and numerical calculations showed that increasing the external 
load on a specimen a transition takes place from a partially failed state of infinite lifetime to a state 
where global failure occurs at a finite time. Two universality classes of creep rupture were identified 
depending on the range of interaction of fibers: in the mean field limit the transition between the two 
states is continuous characterized by power law divergences, while for local interactions it becomes 
abrupt with no scaling. Varying the range of interaction a sharp transition is revealed between the 
mean field and short range regimes. The creeping system evolves into a macroscopic stationary state 
accompanied by the emergence of a power law distribution of inter-event times of the microscopic 
relaxation process, which indicates self organized criticality in creep. 
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Time evolution of physical systems under a steady ex- 
ternal driving is abundant in nature. Examples can be 
found in seemingly diverse systems such as the behav- 
ior of domain walls in magnets [1], earthquake dynam- 
ics [2,3] and creep rupture of solids [4-9]. Recent inten- 
sive experimental and theoretical studies revealed that, 
in spite of the diversity of the governing physical dynam- 
ics, the behavior of these systems share several universal 
features: local stresses arising due to the external driv- 
ing are resolved in the form of avalanches of microscopic 
relaxation events. As a result, the system attains a sta- 
tionary macroscopic state characterized by a scale invari- 
ant microscopic activity. In particular, large efforts have 
been devoted to the study of material failure occurring 
under various loading conditions [4-9] since these stud- 
ies also provide direct information on the dynamics of 
earthquakes [2,3,8,9]. 

Under high steady stresses, materials may undergo 
time dependent deformation resulting in failure called 
creep rupture which limits their lifetime, and hence, has 
a high impact on their applicability in construction ele- 
ments. Creep failure tests are usually performed under 
uniaxial tensile loading when the specimen is subjected 
to a constant load a and the time evolution of the dam- 
age process is followed by recording the strain e of the 
specimen and the acoustic signals emitted by microscopic 
failure events. In spite of the large amount of experi- 
mental results accumulated a comprehensive theoretical 
picture of creep rupture is still lacking. 

In this Letter we study the creep rupture of materi- 
als by means of a novel fiber bundle model. During the 



past years Fiber Bundle Models (FBMs) [10-17] of mate- 
rials damage played a very important role not only in the 
study of fracture but they turned out to be also one of 
the most promising approaches to earthquake predictions 
[3,18]. Our model realistically describes the interaction of 
fibers covering all cases relevant to real materials, further- 
more, it combines their viscoelastic constitutive behav- 
ior and breaking. Analytical and numerical calculations 
show that there exists a critical load that determines the 
final state of the material. Strikingly, we find that in 
creep rupture there are only two universality classes and 
that the distribution of times between breaking events is 
reminiscent of the Omori's law in earthquake dynamics 
[19]. The results suggest that the same dynamics holds 
from the scale of earthquakes down to the rheological 
response of solids. The relevance of our results to exper- 
iments is discussed. 

Our model consists of N parallel fibers arranged on a 
square lattice of side length L having viscoelastic consti- 
tutive behavior [13]. For simplicity, the pure viscoelastic 
behavior of fibers is modeled by a Kelvin- Voigt element 
which consists of a spring and a dashpot in parallel and 
results in the constitutive equation a a = fie + Es, where 
tr is the imposed load, (3 denotes the damping coeffi- 
cient, and E the Young modulus of fibers, respectively. 
In order to capture failure in the model a strain con- 
trolled breaking criterion is imposed, i.e. a fiber fails 
during the time evolution of the system when its strain 
exceeds a breaking threshold e$ , i — 1, . . . , N drawn 
from a probability distribution P(e) = J* Q e p(x)dx. When 
a fiber fails its load is redistributed to the intact fibers 
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according to the interaction law of fibers. It was shown 
in Ref. [13] that in a viscoelastic bundle the fibers break 
one-by-one which implies that such a system is very sen- 
sitive to the details of load sharing, much more than a 
static fiber bundle [10-16] where a large number of fibers 
can break at once in the form of bursts [10]. In order 
to realistically model the stress transfer between fibers, 
recently an adjustable stress transfer function was intro- 
duced, which interpolates between the two limiting cases 
of global and local load sharing [11,12] (GLS and LLS). 
Motivated by results of fracture mechanics we assume 
that the additional load a a dd received by an intact fiber i 
on a square lattice after the failure of fiber j depends on 
their distance and has the form <r a dd — Zr^p . Here 
the normalization factor Z is simply Z = ^2 ieI and 
the sum runs over the set / of intact fibers [15]. The ex- 
ponent 7 is an adjustable parameter of the model, which 
controls the effective range of load redistribution. The 
limiting cases 7 = and 7 — > 00 recover the global and 
local load sharing widely studied in the literature, fur- 
thermore, intermediate values of 7 interpolate between 
them. A comprehensive study of the quasistatic fracture 
of fiber bundles in terms of the adjustable load sharing 
function has been presented in Ref. [15]. 

In the mean field limit, i.e. global load sharing ob- 
tained for 7 = 0, many of the quantities describing the 
behavior of the system can be obtained analytically. In 
this case the time evolution of the system under a steady 
external load a a is described by the differential equation 



where the viscoelastic behavior of fibers is coupled to the 
failure of fibers [13]. For the behavior of the solutions 
e(t) of Eq. (1) two distinct regimes can be distinguished 
depending on the value of the external load <r : When 
a falls below a critical value <r c Eq. (1) has a station- 
ary solution e s , which can be obtained by setting i = 0, 
i.e. a a = Ee s [l — P(e s )\. It means that until this equa- 
tion can be solved for e s at a given external load a Q , the 
solution e(t) of Eq. (1) converges to e s when t —> 00, 
and the system suffers only a partial failure. However, 
when cr exceeds the critical value a c no stationary solu- 
tion exists, furthermore, e remains always positive, which 
implies that for a > a c the strain of the system s(t) 
monotonically increases until the system fails globally at 
a finite time tf. The behavior of e{t) is illustrated in Fig. 
1 for two values of a Q below and above a c . It follows from 
the above argument that the critical value of the load a c 
is the static fracture strength of the bundle [11]. 

The creep rupture of the viscoelastic bundle can be 
interpreted so that for cr < cr c the bundle is partially 
damaged implying an infinite lifetime tf = 00 and the 
emergence of a stationary macroscopic state, while above 
the critical load a a > a c global failure occurs at a finite 
time tf, but in the vicinity of a c the global failure is pre- 
ceded by a long lived stationary state. The nature of the 



transition occurring at a c can be characterized by analyz- 
ing how the creeping system behaves when approaching 
the critical load both from below and above. 

For cr < cr c the fiber bundle relaxes to the stationary 
deformation e s through a gradually decreasing breaking 
activity. It can be shown analytically that e(t) has an 
exponential relaxation to e s with a characteristic time 
scale r that depends on the external load (Tq as 

]/2 

T ~ (<7 C - (T ) , for (T < cr c , (2) 

i.e., when approaching the critical point from below the 
characteristic time of the relaxation to the stationary 
state diverges according to a universal power law with an 
exponent —1/2 independent of the form of disorder dis- 
tribution P [20]. Above the critical point the lifetime tf 
defines the characteristic time scale of the system which 
can be cast in the form [13] 

*/ ~ ( a o - c c )~ 1/2 , for cr > cr c , (3) 

so that tf also has a power law divergence at a c with 
a universal exponent — | like r below the critical point. 
Hence, for global load sharing 7 = 0, the system exhibits 
scaling behavior on both sides of the critical point indi- 
cating a continuous transition at the critical load a c . 

In a creeping system, due to the steady external driv- 
ing, local overloads build up slowly on the microscopic 
level when the breaking threshold of fibers is exceeded 
by the local deformation. The sudden breaking of fibers 
occurring on a time scale much shorter than the scale of 
the driving provides the relaxation mechanism which re- 
solves the overloads in the system. This mechanism con- 
sists of sequential fiber breakings that form an avalanche 
which either stops (below the critical point) or continues 
until the whole system gets destroyed (above the critical 
point). 

Fibers fail one-by-one, furthermore, under GLS condi- 
tions, breakings occur in the order of increasing breaking 
thresholds Si and the time Ai(£j,e i+ i) elapsed between 
the breaking of i-th and i + 1-th fibers can be analyti- 
cally obtained. The inter-event time At is a fluctuating 
quantity which depends both on the breaking thresholds 
and the load level, as illustrated in the inset of Fig. 1. It 
can be observed in the figure that before and after the 
plateau of s(t) the inter-event times are relatively short 
(low peaks), while along the plateau At is scattered over 
a broad interval. The statistics of inter-event times char- 
acterized by the distribution f(At) provides information 
on the microscopic dynamics of creep. f(At) is presented 
in Fig. 2 for a system of N — 2 x 10 7 fibers. Simulations 
revealed that /(At) exhibits a power law of the form 
f(At) ~ At~ b both below and above the critical point 
whenever the macroscopic stationary state characterized 
by the plateau of e{t) is attained. However, the value of 
the exponent b is different on the two sides of the critical 
point, i.e. below a c the exponent of the distribution is 
b = 1.95 ± 0.05 independent of a , while above a c we 
obtained b = 1.5 ±0.05. Increasing the load above a c the 
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stationary state gradually disappears implying that the 
power law regime of f(At) preceding the exponential cut- 
off is getting shorter but the exponent remains the same 
(see Fig. 2). It follows that the creeping bundle self- 
organizes into a critical state: the threshold dynamics of 
the system characterized by a separation of time scales 
of the external driving and the relaxation process leads 
to the emergence of a macroscopic stationary state ac- 
companied by power law distributed microscopic events. 

To explore the effect of the details of load redistribu- 
tion on the creep rupture process we studied how the 
behavior of the system changes in the vicinity of the crit- 
ical point when the load sharing gets localized. Simu- 
lations have been performed varying the effective range 
of interaction of fibers by controlling the exponent 7 of 
the load sharing function. The inset of Fig. 3 presents 
the lifetime tf of a bundle of fibers arranged on a square 
lattice of side length L = 101 as a function of the dis- 
tance from the critical point Act = a a — a c for several 
values of the exponent 7. It can be observed that the 
t/(Aer;7) curves form two groups of different functional 
form: The upper group is obtained for < 7 < 1.95 
when the load sharing is global [15] and Eq. (3) holds [13]. 
However, in the lower group, obtained for 7 > 2.9 when 
the load sharing gets localized [15], tf(Aa;j) rapidly 
takes a constant value showing an abrupt transition at 
the critical load <r c with no scaling, reminiscent of a 
first order transition. The results imply the existence 
of two universality classes in creep rupture character- 
ized by a completely global (GLS), or a completely local 
(LLS) behavior depending on the effective range of inter- 
action 7 with a rather sharp transition between them. 
In order to quantify the behavior of tf(Aa;-f) under 
the variation of 7 we calculated the normalized quantity 
S(l) = [*/(7)-*/(7 = 10)]/[t/(7 = 0)-t/(7=10)], 
where t / (7) denotes the value of tf at the smallest value 
of Act used to calculate tf(A<r;j) at a given 7. Fig. 3 
shows that S(j) provides a quantitative description of the 
creep rupture transition in terms of the effective range of 
interaction so that 5(7) takes value unity for the GLS, 
and it has a value close to zero for the LLS class, respec- 
tively. It can also be observed in Fig. 3 that the tran- 
sition between the two universality classes gets sharper 
around j c w 2 with increasing system size. Real materi- 
als described by a finite value of 7 must fall into one of 
the above universality classes. The existence of only two 
universality classes implies that the mean field analytical 
results can be extended beyond 7 = 0, i.e., they apply 
for a wider interaction range which is relevant for real 
materials. 

Extensive simulations revealed that whenever a macro- 
scopic stationary state is attained by the system, the 
distribution of inter-event times follows a power law ir- 
respective of the range of interaction 7. Below the crit- 
ical point the exponent b of the distribution has a value 
b = 1.95 ± 0.05 independent of 7, while above the criti- 
cal point b is different in the two universality classes as 
illustrated by the inset of Fig. 2. In the LLS class the ex- 



ponent b has practically the same value below and above 

The breaking process of fibers occurring in a solid 
under various loading conditions can be monitored by 
acoustic emission techniques which has also been applied 
to study creep rupture. The statistics of inter-event times 
has been studied in various types of materials like wood, 
plaster, basalt, and fiber glass. It was found experimen- 
tally that the distribution of inter-event times always ex- 
hibits a power law behavior, however, the values of b was 
found to depend on the material falling between 1.2 and 
1.9 [4-9] for o > a c . Hence, our theoretical findings are 
in quite reasonable agreement with the available experi- 
mental results. Moreover, the different values of b below 
and above a c predicted by our model for long range inter- 
actions would correspond to different Omori's exponents 
for foreshocks and aftershocks in earthquake dynamics, 
which has also been observed recently [21]. 

In conclusion, we have identified two universality 
classes of creep rupture depending on the range of in- 
teraction of fibers. The critical behavior of the mi- 
crofracturing process can be seen as the result of the self- 
organization of the system into a macroscopic stationary 
state whose duration depends on the external perturba- 
tion (load). In this sense the system is at some point of 
marginal stability jumping form one metastable state to 
another with power-law distributed waiting times. This 
suggests the existence of a critical dynamics underlying 
the process that seems to indicate self-organized criti- 
cality [4,7,22] in creep. Our theoretical results provide a 
consistent explanation of recent experimental findings on 
the damage process of creep rupture. 

Y. M. thanks A. Vespignani for enlighting discussions 
and A. F. Pacheco for useful comments. This work was 
supported by the project SFB381, by the NATO grant 
PST.CLG.977311. F. K acknowledges the kind hospital- 
ity of the Condensed Matter Group of ICTP Trieste, and 
financial support of the Bolyai Janos Fellowship of the 
Hungarian Academy of Sciences and of the Research Con- 
tracts FKFP 0118/2001 and T037212. Y. M. acknowl- 
edges financial support from the Ministerio de Education 
Cultura y Deportes (Spain) and of the Spanish DGICYT 
Project BFM2002-01798. 



[1] S. Zapperi, P. Cizeau, G. Durin, and H. E. Stanley, Phys. 
Rev. B 58, 6353 (1998). 

[2] Statistical Physics of Fracture and Breakdown in Dis- 
ordered Systems. B. K. Chakrabarti, L. G. Benguigui, 
Clarendon Press, Oxford (1997), and references therein. 

[3] D. L. Turcotte,FractaJs and Chaos in Geology and Geo- 
physics. (2nd Ed. Cambridge University Press, Cam- 
bridge, New York, 1997). 



3 



[4] C. Maes, A. Van Moffaert, H. Frederix, and H. Strauven, 

Phys. Rev. B 57, 4987 (1998). 
[5] A. Guarino, A. Garcimartm, and S. Ciliberto, Eur. Phys. 

J. B 6, 13 (1998). 
[6] A. Guarino, S. Ciliberto, A. Garcimartm, M. Zei, and R. 

Scorretti, Eur. Phys. J. B 26, 141 (2002). 
[7] A. Petri, G. Paparo, A. Vespignani, A. Alippi, and M. 

Constantini, Phys. Rev. Lett. 73, 3423 (1994). 
[8] A. Vespignani, A. Petri, A. Alippi, and G. Paparo, Frac- 
tals 3, 839 (1995). 
[9] T. Hirata, J. Geophys. Res. 92, 6215 (1987). 

[10] M. Kloster, A. Hansen, and P. C. Hemmer, Phys. Rev. 
E 56, 2615 (1997). 

[11] J. V. Andersen, D. Sornette, and K.-T. Leung, Phys. 
Rev. Lett. 78, 2140 (1997). 

[12] F. Kun, S. Zapperi, and H. J. Herrmann, European Phys- 
ical Journal B17, 269 (2000). 

[13] R. C. Hidalgo, F. Kun, and H. J. Herrmann, Phys. Rev. 
E 65, 032502 (2002). 

[14] Y. Moreno, J. B. Gomez, A. F. Pacheco, Phys. Rev. Lett. 
85, 2865 (2000). 

[15] R. C. Hidalgo, Y. Moreno, F. Kun, and H. J. Herrmann, 
Phys. Rev. E 65, 046148 (2002). 

[16] S. Pradhan and B. K. Chakrabarti, Phys. Rev. E 65, 
016113 (2001); S. Pradhan, P. Bhattacharyya, and B. K. 
Chakrabarti, Phys. Rev. E 66, 016116 (2002). 

[17] L. Moral, Y. Moreno, J. B. Gomez, A. F. Pacheco, Phys. 
Rev. E 63, 066106 (2001). 

[18] W. I. Newman, A. M. Gabrielov, T. A. Durand, S. L. 
Phoenix, and D. L. Turcotte, Physica D 77, 200 (1994); 
Y. Moreno, A. M. Correig, J. B. Gomez, A. F. Pacheco, 
J. Geophys. Res. B 106, 6609 (2001). 

[19] F. Omori, Rep. Eart. Inv. Comm. 2, 103 (1894); T. Utsu, 
Fac. Science, Hukkaido University 3, 129 (1969). 

[20] Note that a similar power law divergence of the number 
of successive relaxation steps of a dry fiber bundle sub- 
jected to a constant external load was pointed out in Ref. 
[16]. 

[21] A. Helmstetter, D. Sornette, and J.-R. Grasso, preprint 

cond- mat/0205499 (2002). 
[22] P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 

59, 381 (1987). 



t 1 1 1 1 1 1 1 1 1 1 1 r 



0.5 - A* 




1 ' 1 ' 1 ' 1 ' 1 1 1 ' 1 ' L 

10 20 B0 40 50 60 70 



FIG. 1. e(t) for two different values of the external load a) 
a a < a c , b) a a > a c . The insets present the inter-event times 
At at the time of their occurrence t. 




4 



FIG. 2. The distribution of inter-event times At for 7 = 0. 
Power law behavior can be observed over 5 orders of magni- 
tude. Inset: the exponent b of f(At) as a function of 7 for 
a > a c . 
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FIG. 3. Inset: the lifetime tf as a function of a — o c for 
different values of the exponent 7 between and 100. S(-y) is 
presented in the main figure for several system sizes L. 
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